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P 1. Introduction 



ABSTRACT 



Context, p Pictoris b is one of the most studied objects nowadays since it was identified with VLT/NaCo as a bona-fide exoplanet 
with a mass of about 9 times that of Jupiter at an orbital separation of 8-9 AU. The link between the planet and the dusty disk is 
unambiguously attested and this system provides an opportunity to study the disk/planet interactions and to constrain formation and 
evolutionary models of gas giant planets. Still,/? Pictoris b had never been confirmed with other telescopes so far 
Aims. We aimed at an independent confirmation using a different instrument. 

Methods. We retrieved archive images from Gemini South obtained with the instrument NICI, which is designed for high contrast 
imaging. The observations combine coronagraphy and angular differential imaging and were obtained at three epochs in Nov. 2008, 
Dec. 2009 and Dec. 2010. 

Results. We report the detection with NICI of the planet fi Pictoris b in Dec. 2010 images at a separation of 404 + 10 mas and 
PA = 212.1 + 0.7°. It is the first time this planet is observed with a telescope different than the VLT. 

Key words. Stars: individual (J3 Pictoris) - Stars: early-type - Techniques: image processing - Techniques: high angular resolution 
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The A-type star fi Pictoris is the subject of many attention 
since the discovery of a planet orbiting around (Lagrange et al. 
2009a). The first observation from the Very Large Telescope 
in 2003 reported a candidate companion aligned with the debris 
disk at a separation of -0.41" to the North-East (PA = 34.4°). 
The candidate planet, undetected in early 2009 (Lagrange et al. 
2009b), was finally recovered at ~0.3" to the South- West at 
the end of the same year (Lagrange et al. 2010) definitely con- 
firming its companionship with the star and semi-major axis 
smaller than 15 AU. Given the assumption on the age of this 
system (~12Myr) and its distance (19.44 pc), the candidate 
planet was attributed a mass of about 9 M Jup based on evolution- 
ary models (BarafFe et al. 2003) and later confirmed with more 
detailed photometric works (Quanz et al. 2010; Bonnefoy et al. 
201 1; Bonnefoy et al. 2013) and upper limits from radial veloc- 
ity (Lagrange et al. 2012). Owing to systematic follow-up ob- 
servations the orbit was characterized by Chauvin et al. (2012) 
who determined a semi-major axis of 8-9 AU with low eccentric- 
ity (<0.17). Lagrange et al. (2012) demonstrated that the planet 
is inclined with respect to the main disk and therefore likely re- 
sponsible for the observed warp of planetesimals (Mouillet et al. 
1997; Heap et al. 2000). Recently, a global analysis of the planet 
near-IR photometry and careful comparison to models were pre- 
sented in Bonnefoy et al. (2013). Moreover, ft Pictoris b, the 
closest planet to its parent star ever imaged so far, may well have 
been formed via core accretion. 
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So far, any observations of J3 Pictoris b reported in the liter- 
ature were obtained with NaCo (Nasmyth Adaptive Optics Sys- 
tem and Near-Infrared Imager and Spectrograph), the near-IR 
AO-assisted camera at the Very Large Telescope (Rousset et al. 
2003; Lenzen et al. 2003). A single tentative observation at Keck 
was presented in Fitzgerald et al. (2009) but was performed at a 
time where the planet was angularly too close to the star for be- 
ing detectable and at relatively high airmass. In this paper, we 
present the first analysis of data collected at Gemini South and 
retrieved from the Gemini archive. 



2. Observations 

The Near Infrared Coronagraphic Imager (Toomey & Ftaclas 
2003, NICI), installed at the 8-m telescope, Gemini South, is 
designed to take advantage of the Spectral Differential Imaging 
(SDI) for exoplanets detection (Racine et al. 1999). The instru- 
ment is based on a near-IR (l-5[im) dual-band imager in which 
two images are formed simultaneously (owing to a beam split- 
ter) on two separate detectors. The two spectral channels are 
composed with the same filters so that it is possible to combine 
either two spectrally adjacent filters (CH4L with A = 1.578[xm 
and CH4S with X - 1.652[im) for methanated-object detection 
or any near IR broad-band filters. The pixel scale is 17.932 
and 17.970 mas per pixel (as measured for each detector in Nov. 
2010 at an accuracy of 0.1%) hence providing a field of view of 
about 18.4". NICI also allows for Angular Differential Imaging 
(Marois et al. 2006, ADI) in order to improve the final contrast. 

As for the coronagraphic part, several semi-transparent Lyot 
masks (radii ranging from 0.22" to 0.90") are used in combi- 
nation with several pupil stops undersized with respect to the 
telescope pupil (outer edge stopped from 80 to 95% that of the 
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Table 1. Log of observations which indicates for each data set the date of observation (col. 2), the mask radius (col. 3), the filters (col. 4), the UT 
at start and at end of the sequence (col. 5), the exposure time of individual frame (col. 6), the amplitude of the parallactic angle variation (col. 7) 
and the total exposure time integrated (col. 8). 



Set 


Date 


Mask 


Filters 


UT start/end 


Exp. time 


Parallactic angle 


Total exp. time 






[arcsec] 






[s] 


amplitude [deg] 


[s] 


1 


2008.11.22 


0.22 


CH4L/CH4S 


07:23:40/08:45:15 


3.04 


68.6 


5411 


2 


2009.12.02 


0.32 


CH4L/CH4S 


05:21:10/06:28:11 


4.18/2.66 


39.7 


3224 


3 


2009.12.03 


0.32 


CH4L/CH4S 


05:58:58/08:51:50 


3.04 


67.5 


8086 


4 


2009.12.04 


0.32 


-/H 


05:46:38/08:48:19 


9.88 


72.1 


7884 


5 


2010.12.25 


0.22 


Ks/H 


04:24:03/05:05:37 


1.14 


28.7 


684 


6 


2010.12.25 


0.22 


Ks/H 


05:15:01/07:20:34 


0.76 


43.1 


3268 


7 


2010.12.26 


0.22 


CH4L/J 


06:18:42/06:51:54 


0.76 


9.7 


1064 


8 


2010.12.26 


0.22 


CH4L/J 


06:59:19/07:40:54 


1.14 


8.5 


1425 



pupil). In the following we will refer to two coronagraphic con- 
figurations that associate the 0.22" and 0.32" masks to the 95% 
pupil stop. The focal masks are transmitting a fraction of the 
light (-1:300 for the 0.32" mask) at the center and have a flat- 
topped Gaussian profile. The transmission is therefore not uni- 
form in the field near the mask and has to be taken into account 
when it comes to the photometry of faint companions. Each 
mask is deposited on a substrate which carries its own imper- 
fections. 

While SDI and ADI were already implemented on large tele- 
scopes, NICI is certainly the first instrument fully designed for 
this specific purpose (Liu et al. 2010) and has an important role 
to play in preparing the observations with the near to come 
planet finders namely SPHERE (Beuzit et al. 2008) and GPI 
(Macintosh et al. 2008). 

We retrieved three epochs of observations for J3 Pictoris from 
the NICI archive released on Nov. 24th, 2012, corresponding 
to programs GS-2008B-SV-1, GS-2009B-Q-500, GS-2010B-Q- 
501. The list of available data is given in Tab. 1. 

Data reduction was performed in a similar way as for 
NaCo data using our own pipeline (Bonnefoy et al. 2013; 
Boccaletti et al. 2012) and we followed the guidelines presented 
in other papers reporting NICI observations (Biller et al. 2010; 
Wahhaj et al. 2011; Nielsen et al. 2012). Here, the two channels 
are considered separately. The images are dark- subtracted and 
flat fielded. We did not correct for distortion since we are inter- 
ested in a very small field around the star (the planet is expected 
at a separation closer than 0.4" at these epochs) and we did not 
attempt to perform SDI processing when the CH4L/CH4S filters 
are in use. Flat fields are usually not available for the dates of 
observation so we used calibrations taken several months apart. 
This prevents to account for inhomogeneities in the substrate of 
the focal masks as long as we noticed a drift of the mask po- 
sition by several pixels between the flat fields and the science 
data. In addition, background was not subtracted out since we 
are lacking of sky observations. 

Conversely to NaCo observations of j3 Pictoris the data are 
not saved in data cubes but summed, hence totaling long ex- 
posures of the order of 40-60s. While this is not critical with 
respect to the field rotation at the separation of the planet it can 
be an issue for selecting the best frames. Since the observations 
were obtained with a semi-transparent focal Lyot mask, the se- 
lection is based on the frame with the minimal amount of inten- 
sity (integrated over all the pixels) indicative of a better center- 



ing or a good AO correction. For a given sequence of frames, 
we measured the total flux and retained the frames for which the 
flux is lower than three times the standard deviation of the total 
flux above the minimal flux. For the data presented hereafter the 
selection retained more than 80% of the frames. 

As for the registration of the coronagraphic images we pro- 
ceeded in two steps. First, we applied a rough estimation of 
the star position with a gaussian fit of the images thresholded 
at typically a few percent of the maximum flux, as we did 
in Boccaletti et al. (2012). Then, we performed a more pre- 
cise determination of the star location with Moffat fitting of the 
central attenuated spot behind the semi-transparent Lyot mask 
which is known to move linearly with the actual star position 
(Lloyd & Sivaramakrishnan 2005). The coronagraph centering 
is relatively good in 2008 and 2009 data but features a notice- 
able drift in 2010 data which amounts to several pixels (the 0.22" 
mask is about 12 pixels in radius for comparison). The paral- 
lactic angle is calculated from the azimuth and declination of the 
telescope pointing contained in the FITS headers. 

Out of mask unsaturated images of the star used to determine 
the PSF shape and intensity are not available in the archive. To 
mimic a PSF we used the central spot behind the mask as pro- 
posed in the NICI documentation 1 . We measured a FWHM of 
3.6 pixels equivalent to 64 mas in the Ks band. However, this 
spot is not rigorously an attenuated copy of the star since the en- 
ergy in the center results from the addition of the coronagraphic 
diffraction and the semi-transparent mask transmission. Hence, 
the photometric information provided by this signal is difficult to 
calibrate as it depends on the mask centering, the AO correction, 
the mask size and the wavelength. 

Then, we processed the data with a set of ADI algorithms: 
cADI, rADI (Marois et al. 2006), LOCI (Lafreniere et al. 2007) 
and also KLIP which makes use of a Principal Component Anal- 
ysis (Soummer et al. 2012). For a brief description of these 
algorithms and relevant control parameters see Lagrange et al. 
(2012); Boccaletti et al. (2012). Several parameters for these al- 
gorithms were tested to achieve the best detection performance 
(separation criteria are indicated in Tab. 2). 

As a result, the data taken in 2008 do not show any hints of ft 
Pictoris b as expected from NaCo observations (Lagrange et al. 
2009b) the planet being angularly too close to the star (about 
0.2" in projected separation according to the orbit determination 
in Chauvin et al. (2012)). For the 2009 set of data (CH4L/CH4S 

1 http ://w w w. gemini . edu/sciops/instruments/nici/imaging 
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Fig. 1. Final images processed with rADI (left) and LOCI (right) for 
the Ks and H filter (set 6). The arrow points to the planet. The central 
numerical mask is 0.3" in radius so larger than the coronagraphic mask. 
The field of view is 3"x3". 



and H), it is yet undetected since the Lyot mask is too large 
(r=0.32") with respect to the planet separation (p = 0.3" in Dec. 

2009) and/or the halo is saturated (set 4). We note that it was 
recovered to the SW of the star in NaCo maskless images for the 
same epoch (Oct. to Dec. 2009) but at longer wavelengths (Lp 
filter, 3.8yum) thanks to a lower brightness ratio (Lagrange et al. 

2010) . Moreover, the planet is too warm for methane to produce 
an absorption (T=1700K, Bonnefoy et al. 2013) so spectral dif- 
ferential imaging using CH4L/CH4S filters would not help to 
detect it. Finally, the planet ft Pictoris b is detected in the Dec. 
2010 images (sets 5 and 6) in the Ks band at SNR= 11-15 (mea- 
sured on a 3 -pixel radius aperture with rADI and LOCI) although 
the individual frames are near the saturation level within 0.35" 
from the mask (for the shortest integration time of 0.76s). The 
planet is also seen marginally in the H band (SNR= 5-6). Fig- 
ure 1 presents some images obtained with the rADI and LOCI 
algorithms for both the Ks and H bands. The results obtained 
with KLIP are not particularly of good quality (not shown here) 
and are actually strongly affected by the saturation of the frames 
near the edges of the mask. 

3. Photometry and Astrometry 

In this section we focus on the analysis of set 6 which provides 
the best contrast performance at the location of the planet. The 
extraction of the planet photometry (and its error bar) is delicate 
in these data for two reasons, unsaturated PSFs are missing and 
the focal Lyot mask produces an attenuation on the planet flux. 

There is no published or reported calibration for the 0.22" 
mask attenuation profile. But, according to the NICI instrument 
description, it is expected to be a scaled-down version of the 
0.32" mask attenuation. We retrieved measurements obtained in 
CH4L for the 0.32" mask transmission as a function of X, Y coor- 
dinates. First, we derived the mask center coordinates from this 
low resolution profile, averaged the left and right parts and nor- 
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Fig. 2. Off-axis transmissions of the mask derived from the 0.32" mask 
profile measured in CH4L (solid line) and estimated for the 0.22"/CH4L 
(dotted) and the Ks filter (dashed). 



malized to the intensity in the 0.5-0.8" region where the mask is 
expected to have no influence. Then, we assumed a proportional 
scaling with respect to the size of the mask (0.22/0.32) as well 
as the wavelength (CH4L/Ks = 2.15/1.652), the mask appear- 
ing like if it were smaller at longer wavelengths. The profiles 
we derived are shown in Fig. 2. Assuming the separation was 
390 ± 10 mas (Chauvin et al. 2012) in Dec. 2010, we measured 
that the radial attenuation reaches 0.80 ± 0.02 at the location of 
(3 Pictoris b. As for the central attenuation of a point source, is 
no longer scalable to other wavelengths or mask sizes because 
of the diffraction by the coronagraphic mask. Instead, to mea- 
sure this central attenuation, we used a binary star observed in 
Ks in and out of the 0.22" mask (same object as in Biller et al. 
(2010)) and found a value of 120 ± 5. We retained this calibra- 
tion to derive the planet photometry and we note that his number 
is very different from the attenuation given in the NICI docu- 
mentation (about 300). In addition, it also depends on the star 
centering and AO correction which are not necessarily identical 
when observing the target star. Such variations make the use 
of the central coronagraphic spot a poor photometric calibrator. 
Overall, the uncertainties on the determination of the mask radial 
transmission and central attenuation contributes at a level of 8- 
10% so 0.14 mag hence comparable to the estimation of 1 1% by 
Biller et al. (2010) for PZ Tel b with the same mask/filter config- 
uration (PZ Tel b being at about the same separation as J3 Pictoris 
b). 

To extract the photometry of the planet we proceeded in the 
same way as in Bonnefoy et al. (2011); Bonnefoy et al. (2013) 
so that it is directly comparable to the values derived with NaCo. 
Fake planets (attenuated copy of the PSF) were injected at the 
same radius as the planet (-0.4") but different PAs so that the 
speckle noise is averaged. We considered seven positive fake 
planets (FPPOS) equally separated by 45° and one negative 
fake planet (FPNEG) superimposed with ft Pictoris b. We mea- 
sured the average intensity ratio between positive fake planets 
and the true planet (RPOS) as well as the residual intensity 
at the true planet location when using a negative fake planet 
(RNEG). As we scanned the fake planet to star brightness ra- 
tio from 2.10" 5 to 2.10" 4 , we search for both the value of RPOS 
closer to unity and the minimal value of RNEG. The method was 
repeated for several rADI and LOCI reductions and the mea- 
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Table 2. Signal to noise ratio, magnitude differences in Ks (set 6) and astrometry ( p, PA), measured at the planet location for several algorithms. 



Algorithm 


SNR 


FPOS 


FPNEG 


P 


PA 






(mag) 


(mag) 


(mas) 


(°) 


rADI 0.50 


13.3 


8.71 + 0.49 


8.32 + 0.58 


417+10 


210.98 + 0.70 


rADI 1.00 


15.6 


8.50 + 0.49 


8.57 + 0.58 


402 + 10 


212.30 + 0.70 


rADI 1.50 


11.3 


8.72 + 0.49 


8.77 + 0.58 


402 + 10 


212.28 + 0.70 


LOCI 0.25 


7.9 


8.82 + 0.49 


8.70 + 0.66 


402 + 10 


212.30 + 0.70 


LOCI 0.50 


9.9 


9.06 + 0.49 


9.24 + 0.66 


402 + 10 


212.31 + 0.70 


LOCI 0.75 


8.2 


8.95 + 0.49 


9.25 + 0.66 


401 + 10 


212.44 + 0.70 



sured brightness ratios are reported in Tab. 2 together with the 
SNRs. As for the determination of uncertainties, in addition to 
the errors introduced by the un-perfect knowledge on the mask 
transmission (0.14 mag), we also took into account the observed 
PSF photometric variations (0.23 mag) as well as the disper- 
sion of our measurements (0.12 mag for rADI/FPPOS, 0.21 mag 
for rADI/FPNEG, 0.1 2 mag for LOCI/FPPOS and 0.29 mag for 
LOCI/FPNEG). To be conservative the errors were added lin- 
early amounting to the final error bars indicated in Tab. 2. 

Therefore, the derived magnitude differences of Tab. 2 are 
only in marginal agreement with respect to those measured with 
NaCo (Am Ks = 9.2 + 0.2). While this is compatible within the 
error bars we attribute the discrepancy to 1 / the lack of an unoc- 
culted, unsaturated PSF for the target star, and 2/ from the satu- 
ration of the coronagraphic halo near the position of the planet. 

As for the astrometry, we assumed the star position is deter- 
mined by the central attenuated spot behind the coronagraphic 
mask and we performed Moffat fitting in a small 5 -pixel radius 
aperture centered at the guessed planet location. The obtained 
values (Tab. 2): p = 404+ lOmas and PA = 212.1+0.7° on aver- 
age, are in perfect agreement with those found by Chauvin et al. 
(2012) in Nov. 2010 (390 + lOmas, PA = 212.34 + 2.13°), so 
one month earlier than the present data. The True North position 
is known at an accuracy better than 0.1° (included in our error 
bar) according to the NICI documentation. 

4. Conclusion 

The NICI data released on Nov. 24th, 2012 allowed us to re- 
cover the planet f3 Pictoris b at 404 + 10 mas from the star in the 
South West direction from Dec. 2010 data. It is the first time J3 
Pictoris b is unambiguously detected with a telescope and an in- 
strument which is different than the discovery instrument NaCo. 
Although the presence of the planet is already well attested from 
NaCo observations owing to careful follow-up, the present re- 
sult provides an healthy independent confirmation and demon- 
strates the importance of archiving high contrast imaging data as 
already shown for the HR8799 planets (Lafreniere et al. 2009; 
Soummer et al. 2011; Currie et al. 2012). The astrometry of the 
companion is in conformity with respect to NaCo observations. 
However, we note that the lack of a true PSF of j3 Pic prevents 
us to compare accurately the photometry with previous estima- 
tions from NaCo. More observations with NICI are desirable 
to collect additional photometric measurements and then further 
explore the atmospheric characterization of this planet. With the 
upcoming next generation of planet finders and the hopefully 
large number of planetary mass companion detection, it will be 
mandatory to perform multi-telescope observations for indepen- 
dent confirmations. 

Acknowledgements. We would like to thank the anonymous referee for provid- 
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